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Abstract  -  In  medical  imaging,  a  3D  object  must  often  be 
reconstructed  from  serial  cross-sections.  Usually,  the  cross- 
sections  are  not  closely  spaced  so  that  interpolation  is  needed  to 
solve  the  problem  of  anisotropy.  Different  interpolation 
technique  exist,  applied  either  on  slices  either  on  objects. 

Recently,  a  new  contour  interpolation  method  based  on  a 
zooming  transform  and  on  mathematical  morphology  has  been 
developed.  Thanks  to  the  zooming  transform,  it  works  in  any 
case.  But,  the  obtained  results  are  bigger  and  not  so  smooth  as 
one  could  wish.  This  paper  proposes  a  correction  of  the  method 
in  order  to  avoid  so  big  results. 
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I.  Introduction 

In  many  medical  applications,  a  3D  object  must  be 
reconstructed  from  serial  cross-sections,  either  to  aid  in  the 
comprehension  of  the  object's  structure  or  to  facilitate  its 
automatic  manipulation  and  analysis.  Major  applications  of 
this  technique  can  be  found  in  the  study  of  the  structure  of 
biological  specimens  [1],  the  study  of  the  structural  and 
morphological  characteristics  of  organs  in  biomedical 
sciences  [2],  radiation  treatment  planning  and  surgical 
planning  [3]. 

Usually,  the  medical  imagery  system  do  not  provide 
closely  spaced  cross-sections  so  that  the  spacing  between  the 
slices  is  much  greater  than  the  size  of  a  pixel  within  the  slice. 
Therefore,  establishing  some  kind  of  interpolation  between 
the  slices  is  of  vital  importance  for  3D  reconstruction. 

This  interpolation  step  can  be  reached  either  by  slice 
interpolation  either  by  object  interpolation.  The  slice 
interpolation  was  developed  at  the  beginning  of  the  80's  [4] 
but  object  interpolation  is  a  more  recent  research  axis  (end  of 
the  80's)  which  different  groups  are  interested  in  [5-7]. 

Recently,  a  new  contour  interpolation  method  based  on  a 
zooming  technique  and  on  mathematical  morphology  has 
been  developed  [8].  Thanks  to  the  zooming  transform,  it  is  a 
very  interesting  method  because  it  works  in  any  case,  even 
with  very  complex  shapes  or  with  non  overlapping  contours 
or  non  crossing  contours.  But  unfortunately,  it  gives  bigger 
and  not  so  smooth  results  as  one  could  wish.  This  paper 
presents  a  correction  of  the  method  in  the  calculation  of  the 
inverse  zooming  factor  in  order  to  avoid  so  big  results. 


It.  Method 


A.  Principle 


Considering  the  source  closed  region  Rs  delimited  by  the 
source  contour  Cg  and  the  target  closed  region  RT  delimited 
by  the  target  contour  CT ,  the  principle  of  the  proposed 
method  is  to  apply  a  geodesic  morphological  dilation  of  the 
source  region  inside  a  zoomed  version  of  the  target  region. 
Then,  interpolated  regions  are  computed  by  an  appropriate 
reduction.  This  way  of  enlarging,  interpolation  processing, 
and  reducing  is  called  the  zooming  technique.  The  proposed 
method  can  be  resumed  in  5  steps  : 

1)  operate  a  centroid  registration  transform  which  consists  in 
superimposing  Rg  and  RT ,  according  their  centers  of 
gravity, 

2)  zooming  of  RT  in  such  a  way  it  includes  Rs , 

3)  process  iteratively  a  geodesic  dilation  of  Rs  (until  it 
corresponds  to  the  zoomed  version  of  RT ), 

4)  inverse  zooming  of  each  result  region  of  the  iterative 
process, 

5)  operate  the  inverse  centroid  registration  transform. 

B.  Centroid  Registration  Transform 


The  aim  of  the  centroid  registration  is  to  take  into  account 
the  object  translation  for  slice  to  slice  in  order  to  compensate 
for  large  shifts  in  object  position.  Firstly,  the  centroid  eg  and 
Cj  of  each  region  Rg  and  Rj  are  computed.  Then,  Rs  and 


Rt  are  transformed  by  translation  in  Rg  and  RT 
respectively  in  order  to  superimpose  the  centroids  cs  and  cT 
on  a  unique  point  O  of  the  image  [8]  : 

(  -O  ( 


Rs-T  Rs>cs° 


Rt  =  T] 


Rx  ,cxO 


where  T 


R,U 

v  J 


— ^ 

vector  U  . 


denotes  the  translation  of  a  region  R  with  a 
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C.  Zooming  of  the  target  region 

The  zooming  transform  consists  in  applying  a  scale 
coefficient  or  zooming  factor  to  the  region  RT  in  order  to 

obtain  the  region  RT  in  such  a  way  that  it  includes  the  region 
Rs .  Let  Ls  (respectively  LT )  be  the  longer  segment 
between  eg  (respectively  cx )  and  a  point  belonging  to  Cg 
(respectively  CT )  [8],  Let  Ss  (respectively  ST)  be  the 
shorter  segment  between  cs  (respectively  cT)  and  a  point 
belonging  to  Cg  (respectively  Cx)  [8].  The  zooming  factor 
K  is  computed  as  follows  : 


ST 


This  process  defines  the  transformation  Z  which  leads  to  the 
computation  of : 

R*t  =  z|rt,k| 

where  Z(R,t)  denotes  the  enlargement  provided  by  the 
zooming  transform  with  a  factor  x. 

D.  Interpolation  process 

The  main  stage  of  the  presented  interpolation  method  is 

the  iterative  process  which  deforms  progressively  Rs  into 

' — '  ♦ 

Rt  .  It  is  raised  by  applying  a  geodesic  morphological 
dilation  of  Rs  inside  RT  . 

Let  Bx(x,r)  be  a  geodesic  neighbourhood  ball  centred  at 

x  and  having  r  for  radius,  considered  in  the  region  X  of  the 
plane  ;  it  is  defined  as  : 

Bx(x>r)  =  {yeX,dx(x,y)<r} 
where  dx(x,  y)  denotes  the  geodesic  distance  between  x  and 

y  in  the  region  X.  Then,  the  geodesic  dilation  of  Y  in  X  with  a 
radius  r  is  defined  as  : 

Dx(Y)  =  {x,Bx(x,r)  n  Y  *  0}  =  {x,dx(x,Y)  <  r} 

~  *  ~ 

In  our  case,  X=RT  ,  Y=RS  andr=l. 

Let  N  R  be  the  number  of  iterations  taken  by  the  iterative 

~  *  ~  * 

process  to  reach  RT,  NR  regions  Rj  i=l,  NR  are 
obtained.  For  a  given  a  (0<a<l),  the  corresponding  region  is 

j|C  ' - -  :js  Sfc 

Ra  and  depends  on  a  ( R0  =  Rs  and  Rj  =  RT).  It  is  given 

by  :  Ra  =  R  i(a)  with  l(a)  =  E(aNR)  where  E(.)  denotes 
the  integer  part  ftmction. 

This  process,  described  above,  defines  the  interpolation  IR 
for  the  computation  of  Ra  from  Rs,Rx,a  : 

Ra  =  ^Il(RS’RT’a) 


E.  Inverse  zooming  of  the  selected  dilated  regions 

After  the  interpolation  step,  each  region  Ra  must  be 

resized  correctly  in  function  of  a.  It  is  the  object  of  the 
inverse  zooming  transform  which  consists  in  a  applying  a 
scale  coefficient  or  inverse  zooming  factor  J  given  by  : 

J(a)=^a+1  (*) 

The  region  Ra  is  then  obtained  from  Ra  by  : 

Ra  =  z|Ror,jj 

F.  Inverse  Centroid  Registration  Transform 

After  the  interpolation  process,  an  interpolated  region  is 
obtained  for  each  a.  Then,  each  real  interpolated  region  Ra 
is  computed  by  taking  into  account  the  inverse  centroid 
registration  transform  which  consists  in  computing  the  right 
position  of  the  interpolated  region  : 

(~  -»  -»  'j 

Ra  =  T^Ra,(l-a)Ocs  +  aOcTJ 

III.  Limitation 

Thanks  to  the  zooming  transform,  the  described  method  is 
interesting  because  it  works  in  any  case  even  with  complex 
shapes.  Moreover,  it  does  not  suffer  from  fixed  points 
phenomenon  or  non  crossing  contours.  But  unfortunately,  the 
obtained  results  are  bigger  and  not  so  smooth  as  one  could 
wish  as  illustrated  in  fig.  1  and  fig.  2.  Fig.  1  represents  two 
original  contours  to  be  interpolated,  superimposed  on  a 
unique  image.  Fig.  2  represents  the  superimposed  results  after 
the  computation  of  three  interpolated  contours. 


Fig.  1  :  two  original  contours  superimposed  on  the  same  image 
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Fig.  2  :  Superimposition  of  the  results  obtained  with  the  original  method 


Fig.  3  :  Example  of  two  identical  cercles 
and  an  interpolated  contour  at  a=l/2 


The  main  reason  is  the  unique  speed  of  deformation  for 
each  pair  of  implicitly  matched  points  along  the  geodesic 
paths  reaching  them.  Thus,  the  target  contour  is  not  reached  at 
each  point  at  the  same  time  instead  of  being  reached  at  each 
point  at  the  last  time  of  the  iterative  process  [8]. 

An  other  reason  which  explains  the  size  problem  is  that 
the  calculation  of  the  inverse  zooming  factor  is  not  satisfying 
and  the  object  of  this  paper  is  to  propose  an  other  formula  to 
compute  it. 


IV.  Correction  and  results 


The  inverse  zooming  factor  must  verify  J(0)=1  and 
J  (1)=1  /  K .  In  the  previous  version,  the  inverse  zooming 
factor  was  searched  as  J(a)=Aa  +  B  and  led  to  formula  (*). 

But,  it  is  incorrect.  The  best  way  to  proof  it  is  to  take  a 
simple  example.  Suppose  two  identical  circles  as  Cg  and 
Cj  .  Suppose  their  radius  equal  to  1  and  that  the  zooming 
factor  K  is  chosen  to  be  2.  If  a  =  1/  2 ,  the  corresponding 
contour  after  the  iterative  process  is  a  circle  with  a  radius 
equal  to  3/2  (fig.  3).  With  the  previous  formula  (*),  the 
inverse  zooming  factor  J  for  a  =  1  /  2  is  3/4  instead  of 
2/3. 

To  propose  an  other  calculation  of  J,  let  us  consider  the 
region  Ra  for  a  given  a.  This  region  can  be  expressed  as 


follows  :  =  (1  -  a)Rs  +  aR^ 


Thus,  the  zooming  factor  K(a)  corresponding  to  Ra  is 
given  by:  K(a)  =  (1-  a)  +  aK 

and  the  inverse  zooming  factor  defined  by  J(a)  =  1  /  K(a)  is 


then  : 


J(a) 


1 

(1  -  a)  +  aK 


Fig.  4  :  Superimposition  of  the  results  obtained  with  the  corrected  method 


It  can  be  verified  that  J(0)=1 ,  J(l)=l  /  K  and  that  for  the 
previous  example  of  the  two  identical  circles  J(1  /  2)=2  /  3 . 

Fig.  4  shows  the  results  obtained  by  taking  into  account 
this  new  formula  (**)  for  the  calculation  of  inverse  zooming 
factor.  The  contours  are  not  so  big  than  with  the  previous 
method. 


V.  Conclusion  and  perspectives 

Recently,  a  method  of  contour  interpolation  based  on  a 
zooming  transform  and  on  geodesic  morphological  dilation 
has  been  developed.  It  gives  interesting  results  because  it 
works  in  any  complex  case  but  the  results  are  bigger  and  not 
so  smooth  as  one  could  wish.  A  correction  is  presented  in  this 
paper  in  order  to  limit  the  drawback  of  big  size  phenomenon. 
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The  interest  of  this  method  is  the  use  of  the  zooming 
transform  which  allows  applying  it  on  any  complex  shape. 
The  perspective  will  consist  in  solving  the  problem  of 
smoothness  by  developing  an  interpolation  process  with 
adapted  speed  of  deformation. 
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